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ABSTRACT : 
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1 .   INTRODUCTION 

The  three-stage  axial  flow  compressor  test  rig  of  the  Turbo- 
Propulsion  Laboratory,  Department  of  Aeronautics,  NPS,  will  be  used  to 
carry  out  a  Navy  sponsored  research  program  on  three-dimensional  flow 
phenomena  in  compressor  bladings.   This  compressor  with  a  tip  diameter 
of  36  in.  and  a  hub/tip  ratio  of  0.6  is  presently  equipped  with  lightly 
loaded  bladings  which  are  of  the  free-vortex  type.   A  50  HP  two-speed 
electric  motor  is  now  directly  coupled  to  the  compressor.   Further 
details  of  the  test  rig  are  given  in  Ref.  1  (See  Bibliography  at  end  of 
text) . 

During  the  first  phase  of  the  research  program,  so-called  symmetrical 
bladings  with  high  aerodynamic  loading  will  be  tested.   Further,  in 
order  to  increase  the  pressure  rise  in  the  compressor,  primarily  to 
improve  the  measuring  accuracy,  the  present  motor  will  be  replaced  by  a 
150  HP  motor  that  was  obtained  from  surplus.   This  motor  has  the  follow- 
ing characteristics: 

Make:   General  Dynamics/Electro  Dynamics,  Bayonne,  N.  J. 

Frame:   505  Y  -  Type  TN 

Serial  No. :   70105116  A6 

Cont.  Load:  150  HP,  1180  rpm,  4A0V/3  phase,  190  Amps 

It  is  planned  also  to  make  tests  with  one,  two,  or  three  stages, 
preferable  in  such  a  manner  that  the  full  power  of  the  motor  can  be 
utilized  in  all  three  cases,  requiring  operations  at  different  speeds. 
It  is  intended  to  make  this  possible  by  arranging  a  V-belt  drive  with 
different  replaceable  pulleys. 

The  objective  of  this  study  is  to  establish  the  aerodynamic  design 
of  the  new  blading,  and  to  determine  the  maximum  speeds  at  which  these 
stages  can  operate  at  a  driving  power  of  150  HP.   These  data  will  then 
be  used  in  a  future  report  to  design  the  V-belt  drive  and  the  nozzle  for 
the  measuring  of  the  compressor  flow  rate,  since  Ref.  1  shows  that  the 
presently  used  measuring  method  gives  doubtful  results. 

The  blade  loadings  will  be  such  that  the  NASA  diffusion  factor  does 
not  exceed  a  value  of  0.5  for  the  rotor,  or  0.6  for  the  stator.   All 
three  stages  will  be  identical.   The  flow  conditions  ahead  of  and  after 
the  blade  rows  will  be  determined  with  the  methods  of  Ref.  2. 


2.   FIRST  APPROXIMATION  OF  THREE-DIMENSIONAL  FLOW  PATTERN. 

For  stages  with  equal  energy  input  from  hub  to  tip  without  imposed 

radial  energy  gradients,  the  change  of  the  axial  velocity  component  in 

radial  direction  can  be  obtained  from  Eq.  16(50)  of  Ref.  2.   For  a  first 

approximation  it  is  assumed  that  the  axial  components  V  1  and  V  0  ahead 

ml      mZ 

of  and  after  the  rotor  blades  are  equal  at  the  mean  radius  R  ,  and  that 

m 

the  radial  shift  of  the  stream  surfaces  can  be  ignored.   The  resulting 

velocity  diagram  of  the  blading  at  the  arithmetic  mean  radius  R  is 

m 

shown  in  Fig.  1.   Figure  2  represents  the  conditions  for  a  cylindrical 
stream  surface  at  an  arbitrary  radius  R  between  hub  and  tip.   The  actual 
velocity  diagram  for  a  stream  surface  that  has  the  radii  R-  and  R_  ahead 
of  and  after  the  rotor, respectively,  is  given  by  Fig.  3.   The  symbols 
used  in  these  figures,  and  in  the  following  derivations,  are  defined  in 
Ref.  2. 

From  Eq.  16  (42)  of  Ref.  2  and  by  Fig.  1 
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From  Eq.  16  (53)  of  Ref.  2 
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From  Eqs.  9,  7,  2,  and  10 
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Thus,  the  axial  velocity  components  V  -  and  V  -  at  different  radius 

ratios  R/R  can  be  determined  with  Eqs.  11  and  12  if  the  angles  B, 
m  im 

and  /30   are  specified  at  the  mean  radius  R  . 
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tan  6. 

wR  -  V  . 

■       ul  * 

R/R 
m 

V  ./wR 
ul   m 

V 
a 

1/uR 
m 

With  Eq.  5 

tan3. 

R     Rm   X 

R     R    2 

■   m 

al/wR 
m 

R 

2R 

m   ■ 

R    +   Rm   X 

R»    IT  (13) 

m     k 


2(Val/*R  > 
m 


tan  a.  ■ 

R      R 

V  ,      R   "  -1   X 
ul      m     R 

V  .      2(V  ,/coR  ) 
al        al   m 

M„(l         aZ       m 

tanD2  (V  ,/u.R .  ) 
al   m 

Similarly,  with  Eq.  6 

WR  -  V  2      R   "  /  X 
tan$2  .  __u_   .  ^___  (14) 

Also  from  Fig.  2  and  Eq.  5 


(15) 


and  with  Eq.  6 
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With  these  flow  angles  the  relative  and  absolute  velocities  W- ,  V  , 

and  V, ,  V.  are  known  as  multiples  of  uR  . 
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3.   PERMISSIBLE  FLOW  DEFLECTIONS 

The  NASA  diffusion  factor  D  of  the  rotor  flow  is  defined  by 
[Ref.  6] 


W        AW 

where  aD  is  the  rotor  solidity.   The  diffusion  factor  D  of  the 
K  s 

stator  flow  is 

V1  AV 
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Assuming  that  the  stator  solidity  a     equals  o  ,  and  that  the  chord 

S  R 

of  the  stator  and  rotor  blades  does  not  change  along  the  radius, 
there  is  „ 
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where  o   is  the  solidity  at  R  . 
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For  the  chosen  conditions   the  diffusion   factors  D_  and  D     at   the  mean 
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radius  R  are  equal,  say,  equal  to  D  ,  where  by  Eqs.  21,  2,  and  4 
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Equation  23  can  be  rearranged  to  calculate  directly  the  flow  angle 

B~   for  specified  values  of  8,  ,  D  ,  and  o  .   Equation  23  rewritten 
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It  is  possible  also  to  determine  the  deflection  angle 
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from  Eq.  23.  This  relation  rearranged 
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For  chosen  values  B,  ,  o    ,  and  D  the  flow  angles  B0  can  therefore 

lm   m      m  zm 

be  determined.   The  angles  B,   and  Bn  then  establish  the  quantities 

lm      Zta 

X  and  Y  of  Eqs.  4  and  11,  respectively.   Introduced  into  Eqs.  10  and 
12  these  values  can  be  used  to  obtain  the  changes  of  the  axial  com- 
ponents V  .  and  V  ,.  in  radial  direction.  With  the  angles  B,  and  B« 
from  Eqs.  13  and  14  it  is  then  possible  to  calculate  the  diffusion 
factors  D   and  D  at  particular  radii  with  Eqs.  21  and  22. 

K        S 

4.   APPROXIMATE  OPERATING  CONDITIONS  OF  COMPRESSOR. 

The  following  relations  will  be  used  to  establish  the  main 

parameters  of  the  blading,  which  later  on  will  be  investigated  in  more 

detail. 

If  it  is  assumed  that  V  is  equal  to  the  average  through-flow 

m 

velocity  in  the  compressor  annulus,  the  volume  flow  rate  Q  is 

Q  -  A  V 
cm 

where  A  is  the  cross-sectional  area  of  the  annulus.  With  Eq.  2 
c 
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c   m  tanB,  +  tanB0 
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The  driving  moment  M  for  one  stage  is 

M  -  mK^aV^-  mK  [ft-lb]  (30) 

Because  of  the  small  pressure  rise  in  the  compressor  the  mass  flow 

o 

rate  m  can  be  taken  as 

m  -  Q  p  (31) 
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where  p   is  the  air  density  at  ambient  conditions.   For  z  stages 

(z  ■  1,  2,  or  3)  the  necessary  horsepower  HP  is 
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With  the  rotational  speed  N  in  rpm 
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For  a  tip  diameter  of  3  feet  and  a  hub/tip  ratio  of  0.6 
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Assuming  that  135  HP  are  absorbed  at  the  design  point,  to  have 
available  150  HP  at  overload  conditions,  the  design  operating  speed 
is  obtained  from 

N  - ^^ x   [rpm]  (35) 

z  (tanBlm  -  tanB2m) ]  3 


If  AP  is  the  rise  of  the  total  pressure  in  the  stage,  and  h 
the  total-to-total  stage  efficiency,  there  is 
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With  Eqs.  30,  31  and  3 
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ForAdriving  power  of  135  HP,  the  pressure  rise  AP  is  obtained 
directly  from  the  second  equality  of  Eq .  36.   With  Eq .  33, 


or 
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5.   DETERMINATION  OF  BLADING  PARAMETERS. 

Most  desirable  for  the  present  tests  are  bladings  that  produce 
the  highest  pressure  rise  at  the  lowest  speed  of  rotation.   Such 
bladings  will  also  have  low  flow  rates.   Figures  4  and  5  show  that  the 
angles  6,   should  be  as  high  as  possible  for  these  conditions. 
Figure  6  indicates  that  the  speed  of  rotation  is  not  changing  greatly 
with  6,   and  that  its  magnitude  is  primarily  a  function  of  the  dif- 
fusion factor  D  ,  which  should  be  as  high  as  possible. 

Figure  7  shows  the  flow  deflections  AB  =  8,  -  B„  for  different 

°  lm   2m 

angles  8,   and  diffusion  factors  D  at  a  solidity  a     ■  1,  which 
lm  mm 

establish  the  quantities  X  and  Y  of  Eqs.  4  and  11,  respectively.   These 

values  introduced  into  Eq.  12  show  that  the  axial  velocity  V  _ 

becomes  imaginary  for  B,   greater  than  a  particular  angle  &.        which  is 

depending  on  D  .   These  limits  of  8,    for  V  0  *  0  are: 
m  lm       az 


D    =  0.30    0.35    0.40    0.45 

ID 

6  *  -  45.7°   47.2°   47.8°   48.4° 
lm 

Hence  the  flow  angles  8.   cannot  exceed  about  45  . 

Table  I  gives  the  flow  conditions  along  the  radius  for  8.   ■  45 

lm 

and  D  «  0.4  at  the  mean  radius  R  ,  as  calculated  by  the  previously 

established  relations.   It  can  be  noted  that  the  diffusion  factor  D 

R 

at  the  rotor  tip  is  in  excess  of  0.56,  hence  unacceptably  high.   For 

8,   -  45  and  D  «  0.35  the  resulting  data  are  listed  in  Table  II. 
lmm  ° 

Although  D  at  the  rotor  tip  is  about  0.5,  the  large  deceleration  of 

K 

the  axial  velocity  components  at  the  outer  radius  from  V  - /ttR  ■  0.4445 

to  V  0/u)R  -  0.2629, together  with  the  negative   flow  deflection  AB  at 

az   m         ' 

this  station^make  this  blading  undesirable  for  the  present  purposes. 

The  angle  8,   is  therefore  taken  as  40  and  Table  III  lists  the  data 

obtained  for  D  -0.4.   For  this  value  the  diffusion  factor  D„  at  the 
m  R 

rotor  tip  is  larger  than  0.5.   Repeating  the  calculations  for  8,   ■  40 

lm 

and  D  *  0.35  gave  the  flow  conditions  of  Table  IV  which  seem  acceptable 
m 

in  all  respects.   These  parameters  are  therefore  chosen  for  a  more 

detailed  investigation  of  the  flow  conditions  by  taking  account  of  the 

radial  shift  of  the  stream  surfaces  from  rotor  inlet  to  rotor  discharge 

with  the  resulting  velocity  diagram  of  Fig.  3. 

6.   METHOD  FOR  BETTER  APPROXIMATION  OF  FLOW  CONDITIONS. 

The  quantities  listed  in  Table  IV  will  be  called  the  first 

approximation  of  the  flow  conditions.   Figure  8  shows  the  velocity 

ratios  V  ./(wR  )  and  V  _/(wR  )  of  this  approximation  as  functions  of  the 
al    m       aZ    ro 

radius  ratio  R/R  .   These  values  were  obtained  by  assuming  that  V  .  and 

V  ,  at  R  are  equal,  and  identical  with  the  quantity  V  of  Eqs.  8  and  9. 
az     m  m 

Moreover  a  fluid  particle  entering  the  rotor  at  a  radius  R.  is  supposed 
to  leave  the  rotor  at  the  same  radius,  or  R_  «  R.  ■  R.   This  condition 
led  to  Eq.  1  which  formulates  that  all  fluid  particles  receive  the  same 
energy  increase  AH  in  the  rotor.   However  if  R..  and  R.  differ,  this 
requirement  must  be  expressed  by 

AH  -  u)(R_  V  0  -  R.  V  ,)  -  Constant  (39) 

2  u2    1  ul 
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If  it  is  assumed  that  the  energy  increase  AH  in  the  stage  is 
that  of  the  first  approximation,  there  is  from  Eqs.  1  and  3 

AH  -  ooK  -  u  R  2X  (40) 

m 

where,   from  Table   IV, 

X  -  0.36578 

Now,  if  the  tangential  components  V  .  along  the  radius  R.  ahead  of  the 
rotor  are  taken  to  be  those  of  the  first  approximation,  namely,  by  Eq.  5 

V  .  R  ,   R. 

ul  g  _m   X.      1^  _1  .      . 

wR         "  R-    2      2   R  W) 

ml  m 

the  distribution  of  V  _  along  the  radius  R-  after  the  rotor  is  from 
Eqs.  39,  40  and  41, 


Hfc)'! 


Vu2   Rm  \ 
^=12      f 


Thus,  to  satisfy  Eq.  39,  the  components  V  -  can  be  determined  only  if 
the  corresponding  radii  R.  and  R_  of  the  stream  surfaces  are  known. 
The  results  of  the  first  approximation  will  be  used  to  establish  an 
approximate  relationship  between  R..  and  R?  by  applying  the  equation  of 
continuity.   In  particular,  for  a  chosen  value  of  R..  the  corresponding 
radius  R_  is  found  from  the  condition  that  the  same  mass  flow  rate 

has  to  occur  between  R.  and  the  hub  radius  R_,  ,  as  between  the  radius 

2.  In 

R«  and  the  hub  radius  Rlh«   The  adopted  procedure  is  explained  in  the 
next  paragraph.   With  R.  known  as  a  function  of  R..  ,  the  distribution  of 
the  components  V  ?  in  radial  direction  are  obtained  from  Eq.  42.   The 
axial  components  V  _  along  R-  can  then  be  determined  by  solving  Eq.  16 
(41)  of  Ref.  2,  or 


3(V   )2        2V,   3(R,V  ) 

a2  u2      2  u2  (L'x) 

3  H«  R«     3R^ 
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if  the  curvatures  of  the  generatrices  of  the  stream  surfaces  and  radial 
entropy  gradients  are  ignored.   The  data  obtained  in  this  manner  from 
Eq.  43  represent  the  so-called  second  approximation.   Applied  to  the 
equation  of  continuity  they  permit  a  check  of  the  relationship  be- 
tween R..  and  R-.   If  this  interdependence  does  not  correspond  to  that 
obtained  from  the  first  approximation,  it  will  be  necessary  to  carry 
out  further  iterations. 

Evidently,  if  curvature  effects  and  entropy  gradients  are  ignored, 
the  distribution  of  V  .  along  R-  for  V   of  Eq.  41  is  not  affected  by 
the  radial  shift  from  R..  and  R. ,  hence  the  flow  conditions  ahead  of 
the  rotor  are  those  given  in  Table  IV. 

The  relations  for  the  second  and,  if  necessary,  further  approxi- 
mations, can  be  simplified  by  introducing  the  following  dimensionless 
quantities: 


al 


al 


u)R 
m 


a2 


a2 


0)R 
m 


(44) 


ul 


ul 

u)R 


m 


u2 


u2 

u)R 


m 


ml 


ml 
u)R 


m 


m2 


m2 

uR 
m 


Then,  Eq.  42  rewritten  is 


u2 


£  ("^ 


(45) 


and  Eq.  43  becomes 


('.2) 


3r, 


S  +   fr" 
L  r2       ^  2 


x     3(X  +  rxZ; 
2     3rl 


Since  X  is  a  constant, 
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3r„ 


r   *1 
1  JT2     '  "  f(r2>     <*6> 


Hence,  if  r..  is  known  as  a  function  of  r_,  the  function  f(r?)  can  be 
established.   Denoting  the  as  yet  unknown  value  of  V  _   at  r.  «  1  by 
V  _  ,  the  integration  of  Eq.  46  gives 


(O  2  "  (Vm2*)2-/  f<r2)dr2  <47> 

1 
* 
The  quantity  V  ~  must  be  obtained  by  iterations  with  the  help 

of  the  equation  of  continuity.   Figure  4  shows  that  the  pressure  rise 

per  stage  does  not  exceed  16  inches  of  water.   In  a  stage  with  50 

percent  reaction  the  pressure  rise  in  the  rotor  is  then  about  8  inches 

of  water,  giving  a  pressure  ratio  of  about  1.02  for  an  inlet  pressure 

of  14.7  psia.   The  ratio  Po/P-,  of  the  mass  density  after  and  ahead  of 

the  rotor  will  therefore  not  exceed  a  value  of  1.014. 

The  volume  flow  rate  Q.  ahead  of  the  rotor  between  the  hub 

radius  R..,  and  an  arbitrary  radius  R1  is 


Ql  "  2"/Val  dRl 


Rlh 
With   the  quantities  of   Eq.    44   let 

Q.    -   2ttu>  R  3     Q  *  (48) 

1  ml 

by  introducing  the  dimensionless  quantity 

*      ri 

Qx  m    /     ri   Vai  dri  (A9> 

Rlh/R» 

The  total  volume  flow  rate  Q,     ahead  of  the  rotor  is 

lmax 

Q.     -  2™  R  3  Q*  (50) 

lmax         m   lmax 
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where  Q.    is  obtained  from  Eq.  49  for  the  upper  limit  r,  ■  R,  /R 
xlmax  vv  1    it  m 

of  the  integral,  where  R  is  the  tip  radius. 

1t 
Similarly j for  the  conditions  after  the  rotor  the  volume  flow 

rate  Q„  between  R„.  and  R_  is 

z  zn      I 

Q2   =  2iru)Rm3  Q2*  (51) 


with  Q2*  -    /  r2  V2*  dr2  (52) 

R2h/Rm 

The  total  volume  flow  rate  Q.    at  station  2  is 

2max 


Q.     -=  27ru)  R  3  Q_*  (53) 

2max         m    2max 

The  quantity  0~    is  obtained  from  Eq.  52  for  the  upper  limit 

*<>  "  Ro^/R  of  the  integral. 

2    2t  m  ° 

The  equality  of  the  mass  flow  rates  at  stations  (1)  and  (2) 
can  be  expressed  by 


P,  Q    ■  P^  Qo    k  ' 
1  xmax   K2  v2max  B 

where  p. ,  p_  are  the  average  mass  densities  at  (1)  and  (2) ,  and  kg' 

is  the  so-called  blockage  factor  that  takes  account  of  the  increased 
regions  of  reduced  velocity  in  the  so-called  wall  boundary  layers 
near  the  hub  and  the  tip  at  station  (2) ,  compared  to  those  at  station 
(1).   Then,  with  Eqs.  50  and  53 

*     m     ^lmax m     \max_  (54) 

2max  k 

("2/pl)kB' 

Experience  shows  that  the  blockage  factor  k.  '  varies  between  0.96  and 

_   _  B     * 

0.97.  Hence,  since  p?/  p  ¥  1.014,  the  quantity  Q2    will  be 

determined  by  Eq.  54  with1a  factor  k  -  (1.014) (0.965)  -  0.98  for  the 

*  B  * 

value  of  Q,    obtained  from  Eq.  49.   This  value  of  0„    is  used  to 
Umax  ^  2max 
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3r„ 


1  JT2     "  "  f<r2»     <46> 


Hence,  if  r..  is  known  as  a  function  of  r_,  the  function  f(r»)  can  be 

^  *    l 

established.   Denoting  the  as  yet  unknown  value  of  V  _   at  r,  *  1  by 

*  al  I 

V  -  »  tne  integration  of  Eq.  46  gives 


(Va2*)  2  "  K2*)2-/  f(r2>dr2  <47> 

1 

* 

The  quantity  V  .  must  be  obtained  by  iterations  with  the  help 
of  the  equation  of  continuity.   Figure  4  shows  that  the  pressure  rise 
per  stage  does  not  exceed  16  inches  of  water.   In  a  stage  with  50 
percent  reaction  the  pressure  rise  in  the  rotor  is  then  about  8  inches 
of  water,  giving  a  pressure  ratio  of  about  1.02  for  an  inlet  pressure 
of  14.7  psia.   The  ratio  P^/P-,  of  the  mass  density  after  and  ahead  of 
the  rotor  will  therefore  not  exceed  a  value  of  1.014. 

The  volume  flow  rate  0.  ahead  of  the  rotor  between  the  hub 
radius  R..,  and  an  arbitrary  radius  R..  is 

Ql  "  2*AVal  dRl 
Rlh 

With  the  quantities  of  Eq.  44  let 

Q.  «  2ttw  R  3  Q  *  (48) 

l        ml 

by  introducing  the  dimensionless  quantity 

q*     -  f    rx  Va^*   drx  (49) 


Rih/Rm 

The  total  volume  flow  rate  Q.     ahead  of  the  rotor  is 

lmax 

Q,  -     2iro>  R  3  Q*  (50) 

lmax  m       lmax 
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where  Q.    is  obtained  from  Eq.  49  for  the  upper  limit  r,  ■  R.  /R 
imax  M  vv  1    it  m 

of  the  integral,  where  R  is  the  tip  radius. 

1t 

Similarly j for  the  conditions  after  the  rotor  the  volume  flow 

rate  Q„  between  R_,  and  R_  is 

Z  In  I 

Q2  =  2™Rm3  Q2*  (51) 


r2 
with  Q2*  -    /  r2  V2*  dr2  (52) 

R2h/Rm 

The  total  volume  flow  rate  Q.    at  station  2  is 

2max 


Q2„,ax     '     2™R„3     <C  (53) 

The  quantity  Q-  is  obtained   from  Eq .    52   for   the  upper  limit 

r0  "  Ro^/R     of   the  integral. 

z    zt  m 

The  equality  of  the  mass  flow  rates  at  stations  (1)  and  (2) 
can  be  expressed  by 


Pl  Qmax  '  p2  Q2max  kh 

where  p.,  p_  are  the  average  mass  densities  at  (1)  and  (2),  and  kg' 

is  the  so-called  blockage  factor  that  takes  account  of  the  increased 
regions  of  reduced  velocity  in  the  so-called  wall  boundary  layers 
near  the  hub  and  the  tip  at  station  (2) ,  compared  to  those  at  station 
(1).   Then,  with  Eqs.  50  and  53 

*        lmax  Umax  ,,.,* 

Q2„ax  "  — *  -k—  (54) 

Experience  shows  that  the  blockage  factor  k.  '  varies  between  0.96  and 

_   _  B     * 

0.97.   Hence,  since  p2/  p  "S*  1.014,  the  quantity  Q2    will  be 

determined  by  Eq.  54  with1a  factor  kn  -  (1.0l4) (0.965)  -  0.98  for  the 

*  B  * 

value  of  Q,    obtained  from  Eq.  49.   This  value  of  0„    is  used  to 
lmax  ^  2max 
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check  whether  the  distribution  of  V  .   along  r_,  obtained  from  Eq.  47 

with  a  chosen  value  of  V  _  meets  the  continuity  requirements.   To 

this  end  the  ratios  V  „,  thus  determined,  are  introduced  into  Eq.  52 

az 

and  integrated  from  Rou /R  to  R_  /R  .   If  the  resulting  quantity 
0-    differs  from  Q-    /k_  of  Eq.  54,  the  process  must  be  repeated 

with  a  new  value  of  V  ,   until  agreement  is  reached.   However,  the 

m2  t 

final  distribution  of  V  „  along  r_  must  also  produce  the  same  radial 
shift  of  the  streamlines  as  was  determined  from  the  first  approxima- 
tion, and  which  was  used  to  establish  the  function  f(r~)  of  Eq .  46 

*  l 

and  V  0     by  Eq.  47.   The  equation  of  continuity  is  again  used  to  obtain 

r_  as  function  of  r.  from  the  new  V  „   -  distribution.   If  this  re- 

z  I  az 

lationship  between  r~  and  r.  differs  from  the  one  obtained  with  the 

data  of  the  first  approximation,  it  must  be  used  to  establish  a  new 

function  f(r~)  of  Eq.  46  which  has  to  be  introduced  into  Eq.  47  to 

l*  * 

calculate  V  _  .   This  process  requires  the  same  iterations  for  V  . 

as  discussed  above.   The  iterative  procedure  must  be  repeated  until  the 

functions  r~  vs.  r.  assumed  for  the  determination  of  f(r~)  and  that 

obtained  from  the  resulting  profile  of  V  .  are  identical. 

7.   RESULTS  OF  FLOW  CALCULATIONS 

*  * 

Table  V  lists  the  ratios  V  ,   and  V  -   of  Table  IV,  and  additional 

al       az 

values  at  other  radius  ratios,  which  were  used  to  calculate  the 

*       * 
quantities  Q.   and  Q_   in  accordance  with  Eqs.  49  and  52,  by  integra- 
tions with  the  trapezoidal  rule.   It  can  be  noted  that  the  initial 
assumption  V  1  ■  V  _  ■  V  does  not  satisfy  continuity  since  0~    is 
(0.9832)Q   ^       whereas  by  Eq.  54  with  kg  =  0.98  it  should  be  (1.0204) 
Q1  or  equal  to  0.404947. 

The  ratios  Q-./0-    ,  and  Q2/Qo    of  Table  v  are  plotted  in 

Fig.  9.   For  the  condition  Q0/Q«,    ■  Q, /Q,    ,  the  radius  ratio  r0 

Z   Zmax    l  lmax  z 

for  a  stream  surface  that  has  the  radius  ratio  r..  ahead  of  the  rotor 
can  be  read  on  the  upper  horizontal  scale  of  the  figure.   The  relation- 
ship r1  vs.  r«  thus  obtained  is  shown  in  Fig.  10,  which  also  gives  a 
curve  Ar  ■  (r.  -  r?)  as  function  of  r«  in  a  larger  scale  drawn  through 
the  data  points.   From  this  curve  were  obtained  the  derivatives 
3(Ar)/3r?  by  graphical  means.   Then 
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3rl 
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9(l 

r2  + 
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1 

+ 

3(Ar) 

3t2 

>r2 

3r2 

(55) 


The  numerical  values  of  these  derivatives  and  of  r1  for  different 
ratios  r?  are  listed  in  Table  VI,  together  with  the  function  f(r-) 
of  Eq.  46  for  X  -  0.36578. 

The  function  f(r-)  is  plotted  in  Fig.  11.   Interpolated  values 
between  the  calculated  points  were  used  to  integrate  f(r-)  with  Simp- 
son's rule,  i.e.  by  placing  a  parabola  through  three  neighboring 
points.   The  results  of  the  integration  in  accordance  with  Eq.  47  are 

listed  in  Table  VI  and  plotted  in  Fig.  11.   Table  VI  also  shows  the 

* 
calculated  values  of  V  «   in  accordance  with  Eq.  47,  for  several  assumed 

*        aZ  * 

values  of  V  .  ,  and  the  corresponding  magnitudes  of  Q-  of  Eq.  52 

obtained  by  integrations  with  the  trapezoidal  rule.   As  pointed  out 

earlier,  the  value  of  0-  at  r-  «  1.25,  which  is  denoted  by  Q_    , 

should  equal  0.404947.   For  V  -  «  0.838  this  condition  is  closely 

m2  ^ 

realized  and  the  corresponding  values  of  V    are  taken  to  satisfy 

* 
the  continuity  requirement.   For  v  2  *  0.838  the  ratios  0-/Q?    are 

shown  in  Table  VI  also. 

It  must  now  be  checked  whether  the  relation  between  r.  and  r- 

from  the  first  approximation  of  Fig.  9,  which  was  used  to  calculate 

the  velocity  distributions  of  Table  VI,  agrees  with  the  function  r«  vs.  r. 

that  is  obtained  from  the  second  approximation.   Since  the  change  of 

*  * 

V    along  r1 ,  hence  of  V  .   along  r1 ,  is  the  same  for  both  approxima- 
tions, this  check  can  be  made  by  plotting  Q^/Q«,    as  function  of  r-  for 
both  sets  of  data.   Figure  12  shows  a  curve  of  0-/0-    vs.  r«  drawn 
through  the  values  of  Fig.  9  which  are  marked  by  crosses.   The  data 

points  given  by  circles  in  Fig.  12  represent  the  values  of  Table  VI 

* 
for  V  -  -  0.838.   It  can  be  noted  that  these  points  lie  very  closely 

on  the  curve  of  the  first  approximation,  hence  additional  iterations 
are  unnecessary,  and  the  data  inside  the  heavily  framed  portion  of 
Table  VI  can  be  considered  to  represent  the  actual  flow  conditions 
after  the  rotor  for  the  set  of  initial  assumptions;  that  is,  if  curva- 
ture effects  and  entropy  gradients  are  ignored.   This  final  distribu- 

* 
tion  of  V  -   is  also  shown  by  the  dash-dotted  curve  in  Fig.  8. 
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Table  VII  summarizes  the  flow  properties  ahead  of  and  after  the 
rotor  in  accordance  with  the  data  of  Table  VI.   Since  the  velocity 
diagram  of  the  blading  is  that  of  Fig.  3,  there  are,  with  Eq .  41 
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With  Eq.  45  and  from  Fig.  3 
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Table  VII  also  lists  the  NASA  diffusion  factors  of  Eqs.  17  and 
18,  where  the  subscripts  o  in  (D_)   and  (D_)   indicate  that  these 

K  O         oO 

factors  are  for  a  solidity  a  m  a     ■  1  at  R/R  *  1,  and  that  the  chord 

in  m 

of  the  rotor  and  stator  blades  does  not  change  in  radial  direction. 
In  actuality  the  diffusion  factors  should  be  calculated  for  the 
conditions  along  the  stream  surfaces  by  taking  account  of  their 
radial  shift.   Then,  by  Ref.  3,  for  the  rotor 

W0     R.W  .  -  R0W  . 
DR     1        Wx  +  (TR(R1  +  R2)Wl  (67) 

A  similar  expression  is  defined  in  Ref.  3  for  stators.   However, 
in  the  present  case  where  the  radial  shift  is  slight,  the  relations 
of  Eqs.  17  and  18  can  be  used  with  very  small  error.   In  fact,  the 
differences  are  so  insignificant  that  the  blade  profiles  in  stator 
and  rotor  will  be  determined  for  cylindrical  stream  surfaces,  with 
the  actual  flow  angles  of  Table  VII,  instead  of  along  the  actual 
stream  surfaces.   This  procedure  avoids  the  difficulties  associated 
with  the  definition  of  the  blade  shapes  for  manufacturing  purposes. 
If  the  profiles  are  considered  to  lie  on  arbitrary  stream  surfaces, 
the  resulting  blade  shapes  must  later  be  intersected  by  cylinders, 
or  by  planes  tangent  to  them,  to  be  able  to  obtain  profile  coordi- 
nates that  can  be  used  and  checked  during  manufacture  with  reasonable 
cost.   The  available  calculating  methods  for  the  profile  shapes  and 
their  angular  orientation  are  not  accurate  enough  to  warrant  the 
difficulties  associated  with  the  afore-mentioned  procedure,  particu- 
larly not  because  the  theoretical,  axisymmetric, stream  surfaces  are 
very  nearly  cylindrical  in  the  present  case.   It  must  also  be  recog- 
nized that  the  actual  stream  surfaces  in  a  blade  row  are  not  axisym- 
metric, hence  particles  that  move  along  the  two  sides  of  a  blade 
travel  along  different  paths  with  different  radii  so  that  the  flow 
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does  not  have  the  same  two-dimensional  character  that  exists  in 
rectilinear  cascade  test  rigs  which  were  used  to  establish  the  profile 
design  data. 

In  addition  to,  or  as  a  replacement  of  the  diffusion  factor,  the 
so-called  equivalent  diffusion  ratio  D   is  frequently  used  to  estab- 
lish  the  blade  loading  limits  and  the  profile  losses,  especially  if 
the  flow  incidence  angles  differ  from  those  for  minimum  losses.   By 
Ref.  4,  for  rotors  at  the  design  incidence  angle,  and  for  flows  with 
varying  through-flow  components  on  cylindrical  stream  surfaces, 
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(68) 


For  stators  at  the  same  conditions, 


V  -  %    i  i-«  +  If  <»%  5t^  I 


Experience  has  shown  that  the  quantities  D   should  not  exceed  values 
of  about  1.7  since  incipient  flow  separation  is  likely  to  occur  for 
D_.  between  2.0  and  2.2.   The  subscripts  o  in  (D_._)   and  (D_„) 

of  Table  VII  indicate  that  o.  ■  o  «  1  at  the  mean  radius  R  and  that 

R    S  m 

the  blades  have  constant  chord  at  all  radii. 

8.   BLADE  PROFILES  AND  STAGGER  ANGLES 

The  thickness  distribution  of  Fig.  13  will  be  used  for  the 
profiles  of  the  stator  and  rotor  blades.   It  corresponds  to  a  C.A  base 
profile  (See  Ref.  5),  except  for  the  portions  near  the  entrance  which 
were  thinned  to  arrange  a  smaller  leading  edge  radius.   The  camber 
line  of  the  profile  is  supposed  to  be  a  circular  arc.   Figure  14  shows 
the  method  that  will  be  adopted  to  determine  the  profile  coordinates 
for  known  camber  angles  a  and  maximum  thickness  ratios  t/c.   The 
orientation  of  the  blade  profiles  in  the  cascade  is  fixed  by  the  stagger 
angle  y   as  indicated  in  Fig.  14. 
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i  The  profile  data  are  determined  with  the  two-dimensional  method 
of  chapter  VI  of  Ref.  6,  without  taking  account  of  the  corrections 
for  three-dimensional  effects  that  are  given  in  chapter  VII  of  Ref.  6. 
It  is  felt  that  these  corrections  are  not  sufficiently  well  supported 
by  experiments.   Moreover,  since  the  blades  in  the  test  compressor 
can  be  set  at  different  angles  the  twist  of  the  profiles  along  the 
blade  height  is  more  important  than  the  stagger  angle  per  se. 
The  camber  angle  <j>  for  a  rotor  profile  is  obtained  from 

<|>  ■  o    o  (70) 

1  -  m  +  n 


where 


*o  '  (Ki>sh  (Vt  (io>10  (71) 


6o  -  <Vsh  <Vt  (6o>10  (72) 

(i  )  «  zero  camber  incidence  angle  for  t/c  »  0.1 

(6  )  _  ■  zero  camber  deviation  angle  for  t/c  »  0.1 

The  recommended  value  of  the  correction  factors  (K. )  .  and  (K.)  ,  is 

i  sn       o  sn 

1.1  for  C.A  profiles.   The  design  incidence  angle  i  for  minimum  loss  is 

i  «  i  +  n<(>  (73) 

o    r 

and  the  deviation  angle  6  for  these  conditions  is 

6*6     +  m<J)  (74) 

o    T 

The  quantities  (i  )   ,  (6  )   ,  m,  and  n  are  functions  of  the  inlet 

0  10    °  10 
flow  angle  B,  and  the  blade  solidity  a.   The  correction  factor  (K  ) 

t 
and  (K.)   depend  on  the  thickness  ratio  t/c  only.   The  blade  stagger 

angle  y  of  Fig.  14  is  equal  to 


Bx  -  i  -  |  (75) 
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The  stator  blade  profiles  are  obtained  with  Eqs.  70  to  74  also  If 
B,  is  replaced  by  cu,  and  AB  by  Act  ■  a?  -  a.  .   The  blade  stagger  angle 


is  then 


o2  -  i  -  f  (76) 


For  particular  conditions  the  method  of  Ref .  6  gives  excessively 
large  negative  incidence  angles.   Although  negative  incidence  angles 
are  beneficial  for  good  stall  margins,  they  may  cause  increased  losses 
at  the  design  point  because  of  separations  on  the  pressure  side  near 
the  leading  edge,  especially  if  the  latter  has  a  small  radius.   If 
the  above-mentioned  situation  occurs,  reduced  negative  incidence  angle 
i'  will  be  assumed.   Reference  6  shows  that  for  these  conditions  the 
corresponding  deviation  angle  6  is  obtained  from 

«'-«  +  (!'-  i)   (||)  (77) 


where  the  function  (d6/di)  depends  on  solidity  and  inlet  flow  angle. 
For  a  stator  blade  profile,  for  instance,  there  is 

Aa-4+i-6«$'+i'-6' 
where  the  first  expression  on  the  right-hand  side  holds  for  the  design 
conditions  obtained  with  the  usual  approach,  whereas  the  second 
holds  for  an  assumed  incidence  angle  i*  from  which  the  necessary  blade 
camber  angle  <J>*  can  be  determined.   With  Eq.  77 


»         »      » 

■  Act  -  i  +  6   -  Aa 


-  i'  +  6  +  <i'  -  i)   (||)    (78) 


For  this  modified  camber,  the  blade  stagger  angle  y     of  a  stator 

is  y*   »  a,  -  i'  -  ±  (79) 

2 


The  axial  distances  between  the  blade  rows  of  the  compressor 
are  2.875  inches.   The  axial  clearance  between  the  rows  must  be  about 
0.5  in.  to  provide  room  for  flow  survey  probes, so  that  the  axial  blade 
widths  must  be  less  than  about  2.4  inches.   Further,  since  each  rotor 
has  30  and  each  stator  32  blades, the  solidites  of  the  rotor  and  the 
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stator  cascades  must  remain  within  certain  limits. 

The  diffusion  factors  D  and  the  equivalent  diffusion  ratios  D 

of  Table  VII  are  measures  for  the  aerodynamic  loading  of  the  blade 
profiles  in  a  cascade.   Excessive  loading  is  associated  with  flow 
separations,  increased  losses,  and  reduced  efficiencies.   At  low 
blade  loadings  the  efficiency  of  the  cascade  is  decreased  also  because 
of  large  boundary  layer  losses.   Optimum  efficiencies  of  a  row  of  blades 
are  obtained  if  at  all  radii  the  blade  profiles  have  the  highest  possible 
loading  and  the  same  margin  of  safety  against  stalling.   The  choice  of 
the  latter  depends  on  the  operating  characteristics  of  the  system  for 
which  the  compressor  is  used. 

Table  VII  shows  that  this  condition  is  not  realized  in  symmetri- 
cal bladings  with  blades  of  constant  chord.   The  diffusion  factor  of 
the  rotor  is  highest  at  the  outer  radius  and  very  much  below  permissible 
limits  at  the  hub  radius.   The  opposite  trend  would  give  a  better  per- 
formance since  it  is  known  that  for  the  same  diffusion  factor  the  losses 
are  higher  at  the  rotor  tip  than  at  the  hub.   The  data  for  the  stator 
in  Table  VII  show  that  the  profiles  are  lightly  loaded  at  the  outer 
radius  and  have  increased  loading  at  the  hub  radius.   For  stators  with 
radial  clearance  gaps  at  the  hub  radius  the  loading  should  be  reversed 
also  to  obtain  optinum  performance, since  the  stator  losses  near  the  hub 
radius  are  larger  than  those  at  other  radii  even  if  the  diffusion  factor 
were  constant  everywhere  along  the  blade. 

Within  the  limitations  imposed  by  the  rotor  and  stator  blade 

numbers  and  the  necessity  of  maintaining  sufficiently  large  axial  blade 

clearances  for  the  arrangement  of  flow  survey  probes,  the  blade  chords 

of  rotor  and  stator  will  be  varied  along  the  radius  to  slightly  improve 

the  distributions  of  the  diffusion  factors  of  Table  VII.   The  chosen 

chords  for  the  rotor  blades  are  given  in  Table  VIII,  those  for  the 

stator  blades  in  Table  IX.   From  Table  VIII  it  can  be  seen  that,  whereas 

for  constant  chord  the  factor  D  for  the  rotor  changed  from  0.2062  to 

K 

0.4376  from  hub  to  tip,  it  now  has  values  of  0.2498  and  0.4274  at  these 
respective  stations.   The  stator  solidities  obtained  with  the  chosen 
chords  are  generally  lower  than  those  for  constant  chord  in  Table  VII, 
which  gave  variations  of  D  from  0.4107  to  0.3326  from  hub  to  tip.   By  Table 
IX  these  respective  values  are  0.4240  and  0.3762  for  the  stator  blades 
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with  varying  chord.   These  data  show  that  with  the  symmetrical  blading 
of  Fig.  3  the  earlier-mentioned  optimum  blade  loadings  cannot  be  achieved 
by  varying  the  blade  chord.   Symmetrical  bladings  produce  favorable 
conditions  as  far  as  the  uniformity  of  the  inlet  flow  velocities  along 
the  radius  is  concerned.   Table  VII  shows  for  instance  that  the  relative 
rotor  inlet  velocity  VL  decreases  from  100  percent  at  the  hub  to  95.5 
percent  at  the  tip.   The  stator  inlet  velocity  V-  decrease  in  radial 
direction  also,  namely,  from  100  percent  at  the  hub  to  81  percent  at 
the  outer  radius.   The  actual  magnitudes  of  W.  and  V_  are  nearly  equal. 
Hence,  in  contrast  to  other  blading  types,  symmetrical  bladings  are 
not  Mach  number  limited  by  the  conditions  at  the  outer  radius.   For  a 
given  maximum  Mach  number,  they  can  operate  at  higher  peripheral  speeds 
than  other  bladings  and  are  therefore  capable  of  producing  a  higher 
pressure  rise  per  stage. 

The  best  design  of  a  blading  would  however  be  one  that  has  uniform 
and  equal  velocities  W.  and  V.  along  the  radius,  which  are  small  with 
respect  to  the  peripheral  speed,  and  where  the  diffusion  factors  for 
all  profiles  are  as  high  as  the  necessary  stall  margins  permit. 
Investigations  whether  such  bladings  are  possible  will  be  carried  out 
with  a  new  computer  program  that  is  based  on  the  calculating  approach  of 
this  report,  but  with  the  additional  refinement  that  the  effect  of  the 
different  losses  along  the  radius  will  be  taken  into  account. 

The  profile  thicknesses  of  Table  VIII  and  IX  were  chosen  for  practi- 
cal and  aerodynamic  reasons.   For  the  rotor  the  blade  thickness  increases 
about  parabolically  from  tip  to  hub  to  obtain  reduced  centrifugal 
stresses,  and  the  thickness  of  the  stator  blades  varies  about  linearly 
to  obtain  small  ratios  t/c  at  the  hub  radius,  and  sufficient  blade 
thickness  at  the  outer  radius  for  the  attachment  of  the  blade  to  its 
root  section  with  the  circular  platform  that  is  shown  in  Fig.  9  of  Ref.l. 

Whereas  the  incidence  angles  i  of  the  rotor  blade  profiles  of  Table 

VIII  seem  reasonable,  excessively  large  negative  incidence  angle  are 

obtained  with  the  design  method  of  Ref .  6  for  the  stator  profiles  of 

Table  IX  at  large  radii.   The  data  inside  the  heavily  framed  portion  of 

Table  IX  are  therefore  replaced  by  those  for  the  assumed  incidence 

» 
angles  i  ,   The  last  lines  in  Tables  VIII  and  IX  give  the  axial  widths 

of  the  rotor  and  stator  blade  profiles.   For  an  axial  spacing  of  2.875  in. 
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of  the  stator  and  rotor  rows,  the  average  axial  clearance  between  the 

rows  Is  2;875-(l    +1   )/2.   Hence  this  axial  clearance  is  about 
axr    axs 

0.47  in.  at  the  hub  radius  and  0.76  in.  at  the  outer  radius.   These 
clearances  are  large  enough  for  blades  with  the  usual  radial  stacking 
line  through  the  centers  of  gravity  of  the  different  profiles.   It  must 
however  be  examined  whether  for  blades  with  swept-back  leading  edges, 
which  are  to  be  tested  during  the  program,  the  blade  chords  of  Tables 
VIII  and  IX  permit  to  arrange  inter-stage  instrumentation. 

The  data  of  Tables  VIII  and  IX  will  be  used  to  establish  the  profile 
coordinates  with  the  relations  of  Fig.  14.   For  each  profile  there  must 
be  determined  its  cross-sectional  area,  the  coordinates  of  the  center  of 
gravity  and  the  angles  of  the  principal  axes  with  respect  to  the  refer- 
ence line  of  Fig.  14,  and  the  minimum  and  maximum  moments  of  inertia. 
These  quantities  are  necessary  to  establish  the  stresses  and  the  vibra- 
tory characteristics  for  particular  three-dimensional  arrangements  of 
the  profiles  in  the  build-up  of  the  blades. 

9.   COMPRESSOR  PERFORMANCE. 

If  the  losses  in  the  inlet  guide  vanes  of  the  compressor  are 
ignored,  the  total  pressure  P   ahead  of  the  first  rotor  at  station  (1) 

is  everywhere  constant.  However  the  static  pressure  p.  at  station  (1) 
changes  along  the  radius,  since  for  incompressible  flows 

"1  '  Ptl-2V12  (80) 

At  station  (2)  after  the  rotor  the  total  and  static  pressures  are 

denoted  by  P  _  and  p  , 

where 

"2  "  Pt2  "  2  V2  <81> 

In  contrast  to  prior  usage  in  this  report  the  region  after  the  stator 
is  called  station  (3) .  At  this  station  exist  the  same  velocities 
and  flow  angles  as  at  station  (1)  but  the  pressures  differ  from  those 
at  (1),  and  the  static  and  total  pressures  p~  and  P  _  are  related  by 

p3  -P^-fV  '  (82) 
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From  Eq.  7(44)  of  Ref.  2  the  conditions  along  a  streamline  in  the 
rotor  between  stations  (1)  and  (2)  can  be  expressed  by 

P2  +  2  W22  '  2  "'  R22  *  Pl  +  2  Wl2  "  2  ^   "  APR  (83) 

where  AP  represents  the  losses  in  the  rotor.   In  Eq.  83  the  differ- 
ences  in  the  geopotential  energies  of  the  fluis  are  ignored.   It  is 
customary  (See  Ref.  6)  to  express  the  rotor  losses  by 

APR  -  VRfWl2  (84) 

In  Ref.  3  the  pressure  loss  coefficient  Y  is  related  to  the  diffusion 

R 

factor  D  and  the  rotor  solidity  oD  by 
K  R 


2o     r  2.91       2.02      3.77  1 

YR  *  "EoTF     °*00A  +  °-0639(Dr  +  °-1>  +  0.228DR    (1-XR) 


where 

R2t  "  R2      r2t  "  r2 

xr"   r    -  r      =  r^r  (86) 

R    R2t   R2h     r2t   r2h 

The  subscripts  t  and  h  refer  to  the  outer  and  hub  radius,  respectively. 

Equation  85  is  an  approximate  mathematical  formulation  of  the  empirical 

loss  correlations  of  Ref.  6  which  attempt  to  account  for  tip  clearance 

effects  with  the  factor  X  ,  which  increases  the  losses  toward  the  rotor 

K 

blade  tip. 

In  Ref.  6  the  loss  in  total  pressure  in  the  stator  blade  row  is  defined 

by 

APS  -  Pt2  "  Pt3  "  YS  2   V22  (87) 


where  by  Ref.  3  the  stator  loss  coefficient  Y  can  be  approximated  by 


(85) 
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2os 

'  «  x_ 

S   cosa- 


2.91       2.02        3.77] 
0.004  +  0.0639  (Dg  +  0.1)  +  0.057DS    (1  -  Xg)         (88) 


The  quantities  a  and  D  are  the  stator  solidity  and  the  stator  dif- 
fusion factor, respectively ,  and 

XS  "  r1  "  Y  (89) 

S     Rit   Rlh 


Due  to  the  ratio  X_  the  loss  coefficients  Y  of  Eq.  88  increase 
toward  the  hub  radius  where  the  stator  blades  have  a  radial  clearance 
gap.   Evidently  Eqs.  84  and  88  do  not  take  proper  account  of  the  tip 
clearance  and  end  losses.   A  discussion  in  Section  13.7  of  Ref.  2 
shows  that  these  losses  must  depend  on  the  blade  loading,  the  ratio  of 
tip  clearance  and  blade  height,  and  the  ratio  of  blade  spacing  and  blade 
height,  in  addition  to  solidity.   The  present  research  program  endeavors 
to  establish  a  better  and  more  physically  correct  evaluation  of  the 
end  losses  and  tries  to  investigate  whether  peculiar  blade  shapes  can 
be  used  to  minimise  their  detrimental  effects  on  stage  performance. 
From  Fig.  3 

2  2      2      2   2  2 

W.    =   (coR.  -  V  .)   +  V  .    -  co  R.   -  2ojR,V  ,  +  V- 
1         1    ul      al         1       1  ul    1 


W/     -      (coR2  -  Vu2)2  +  Va22     *     u,2R22  -  2WR2Vu2  +  V,,2 

Introduced   in  Eq.  83,    and  with  Eq.    84, 

2  2                                                               2 

P2  +  f  2        -  Pl  +  fx     +  P»(Vu2-  RlVul>    "  YR  fl 


With  Eqs.    39,    40,   80,    and   81 


Pt2-Pu-     pM2Rm2     X-YR|Wl2  (90) 
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where  X  is  a  constant,  equal  to  0.36578, for  the  bladings  under  considera- 
tion.  With  Eq.  64  the  above  relations,  rewritten,  gives 

Pt2  *  Ptl  YR    *  2 

*t2  -  —2 — r   ■  x  -  r  (wi  >  «» 

pu>   R 
m 

From  Eq.  87,  with  Eq.  65 


P    —   P         Y        o 

1£2 lt3   .  _£     *  2  (92) 

2    2         2   v  2  ' 

pw   R 
m 

and 

Pt3  "  Ptl  \  *  2         ts         *  2 

"t3  '    2  .2      "   *  --<«!>-/  <V2  >  (93) 

po)    R 
m 


For  frictionless  conditions  the  coefficients  Y  and  Y_  would  be  zero 

R      S  2  2 

and  the  theoretical  total  pressure  rise  in  the  stage  would  be  Xpw  Rm 
by  Eq.  93.   Hence  the  total-to-total  stage  efficiency  becomes 

,     VV   ^TS(V2>  ,   N 

n  -  1  -  (9A) 

t  2X 

Equation  94  should  actually  be  evaluated  along  particular  stream 

■k  it 

surfaces;  that  is,  W.   and  V_   should  be  taken  at  the  corresponding 
radii  R..  and  R_  of  a  stream  surface.   However,  because  of  the  small 
radial  shift  of  the  flow  paths  in  the  stage,  and  for  the  other  reasons 
explained  in  paragraph  8,  Eq.  94  will  be  determined  for  cylindrical 
stream  surfaces  and  constant  radii  R  =  R.  «  R_.   A  volume  flow  averaged 
stage  efficiency  TP  can  be  evaluated  from  the  different  efficiencies 
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n.    along  the  radius  by 


ru 


/    d"i  "t 


n   .  lih (95) 

rlt 


/ 


^1 

rlh 


With  Eqs.  49  and  50,  with  r,.  -  R,u /R  ,  r_  *  R.^/R  , 

In    In  m   It    It  m 


/ 

J       rl  Val  nt  drl 

"<  -^ s (96) 

lmax 

where  from  Table  V,  Q      -  0.397007.   The  average  total  pressure 
P  _  at  station  3  Is  then  obtained  from 

ft3  "  ?1  ".  2       "  \X  (97) 

pw   R 
m 

In  addition  to  the  changes  of  the  total  pressures  P  _  and  P  _  in 
radial  direction,  the  distributions  of  the  static  pressures  p1 ,  p~, 
and  p_  are  of  interest  also  for  comparisons  with  the  measured 
pressures  and  for  the  determination  of  the  aerodynamic  blade  forces. 
From  Eqs.  80  and  63 

ni-  -V-f   -  -i<V  (98) 

pa>  R 
m 
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From  Eqs.  81,  91,  and  65 


P2  "  Ptl         YR     *  2   1     *  2 

n2 "       2 — r   '    x  -  -r    (W1  )    "  1  (V?  >         (99> 

n 


From  Eqs.  82,  93,  and  63 


P3  -  Ptl        YR    *  2   Y<?   *  2   1    *  2 
n3  -  -^ — r  »  X  -  ^  (Wx  )  -  /(v2  )  -  i(v1  )     (100) 

pu>  R 
m 


Table  X  shows  the  loss  coefficients  in  rotor  and  stator  for  the 
solidities  of  Tables  VIII  and  IX,  and  gives  the  dimensionless  pressure 
coefficients.   It  can  be  noted  that  highest  efficiency  n  -  0.9422 
occurs  at  the  mean  radius,  whereas  the  efficiencies  at  the  tip  and  the 
hub  are  0.8238  and  0.8857,  respectively.   The  last  line  of  Table  X 
lists  the  integrands  of  Eq.  96.   By  trapezoidal  integration  the  inte- 
gral of  Eq.  96  is  0.363035,  hence  the  average  stage  efficiency  is 

0.363035    n  Q  ,, 
nt  "  0.397007  "  u-yi4A 

Then,  by  Eq.  97,  with  X  -  0.36578 


Pt3  -  Ptl 

Ht3  "    2  .  2    "  °'33AA8 
pu   R 
m 


Figure  15  is  a  plot  of  the  pressure  coefficients  of  Table  X  to 

show  the  pressure  changes  along  the  radius  of  the  stage.   It  can  be 

noted  that  equal  increases  of  the  static  pressures  in  rotor  and  stator 

occur  at  a  radius  ratio  r  »  R/R  of  about  1.1.   At  smaller  radii  most 

m 

of  the  static  pressure  rise  occurs  in  the  stator,  whereas  the  opposite 
occurs  between  r  ■  1.1  and  the  tip.   Hence  because  of  the  resulting 
radial  velocity  distributions  a  symmetric  blading  does  not  have  a 
constant  degree  of  reaction  of  50  percent. 
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With  the  established  data  the  operating  conditions  of  the  compressor 
can  be  determined  so  that  135  HP  are  absorbed  by  one,  two,  or  three 
stages.   By  Eqs.  30,  31,  40,  and  50,  the  driving  moment  M  for  one  stage 
is 

M  -  0.     p   wR  2  X  -  2ir  w2R   p   Q,*   X 
lmax  o   m  mo  lmax 

* 
where  0.     -  0.397007  from  Table  V.   Then,  by  Eq.  36,  for  135  HP, 
lmax  _         .. 

R  -  1.2  ft,  and  p   -  2.371(10~J)  slug/ft  , 
m  o 


u3  .  (135X550)   .  13_793  (106, 

2ir  p   R   Q.      X 
o  m   1  max 

or 


en  ■  239.82  radians/sec 


30 
and       N  -  —  u>  *  2290  rpm 


Then,  by  Eq.  50, 


Q.    -  1034   ft3/sec 
xlmax 

From  Eq.  97  with  X  -  0.36578  and  n"t  -  0.9114 

?  0  -  Pt.  -  p  w2R  2  X  IT  -  62.4  lb/ft2  -  12.6  in.  H.O 
tJ    tl    o    m      t  2. 

The  following  table  gives  these  data  if  one,  two,  or  three  stages 
absorb  a  driving  power  of  135  HP. 
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One  Stage  Two  Stages  Three  Stages 
Z-l    Z-2      Z-3 

Speed  N  (rpm)      2290    1818       1588 

Volume  Flow  Rate  Qlmax   (ft3/sec)      1034    820        718 

Pressure  Rise  Z<?t3-*tl)    (in.  H20)      12.6    15.9       18.1 
Rotor  Tip  Speed  uRt       (ft/sec)      359.7   285.5      249. 4 

Before  deciding  on  a  drive  for  all  three  speeds,  a  study  will  be 
undertaken  to  determine  the  blade  and  rotor  stresses,  and  to  investi- 
gate whether  critical  shaft  speeds  and  bearing  capacities  permit 
trouble-free  operation  for  all  cases. 
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FIG.    1     VELOCITY  DIAGRAM  OF  SYMMETRICAL 
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FIG.    2     VELOCITY  DIAGRAM  OF  SYMMETRICAL 
STAGE  FOR  FLOW  ON  CYLINDRICAL 
STREAM  SURFACE  WITH  R?  -  R 
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FIG.  3  VELOCITY  DIAGRAM  OF  STAGE  WITH 
FLOW  ON  ARBITRARY  AXISYMMETRIC 
STREAM  SURFACE 
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FIG.   4     THEORETICAL  PRESSURE  RISE  APt/"H  IN  COMPRESSOR   FOR  POWER  OF   135  HP 
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FIG.    5     VOLUME   FLOW  RATE  Q  OF  COMPRESSOR   FOR  POWER  OF   135  HP 
SYMMETRICAL  BLADING  WITH  a     -   I 
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FIG. 6     SPEED  OF  ROTATION  N  OF  COMPRESSOR   FOR   POWER  OF   133  HP 
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FIG.    7      FLOW  DEFLECTIONS   IN   ROTOR  CASCADE  AT  DIFFERENT  VALUES 

OF  FLOW  INLET  ANGLE        8  AND  NASA  DIFFUSION    FACTOR  D 
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FIG. 11      FUNCTION    fu?)    AND  yf(r   )   dr2   FOR   SECOND  APPROXIMATION 
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FIG.  14     DETERMINATION  OP  PROFILE  COORDINATES 
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FIG. 15  CALCULATED  PRESSURE  DISTRIBUTIONS  IN  STAGE  (SEE  TABLE  X) 

Pt  •  Total  Pressure  ,  p  -  Static  Pressure 

R   -  Radius,  u>  -  Angular  Velocity,  p  -  Mass  Density 

Subscripts:  (*)  Ahead  of  Rotor 

(2)  After  Rotor, Ahead  of  Stator 

(3)  After  Stator 

(m)   Arithmetic  Mean  Radius  of  Blading 
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TABLE    I 

THREE-DIMENSIONAL    FLOW  CONDITIONS    IN    COMPRESSOR    STAGE  WITH  VELOCITY  DIAGRAM 
OF  Fib.  2   ON    CYLINDRICAL   STREAM   SURFACES 


i  AT  1 
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TABLE      I 
THREE-DIMENSIONAL    FLOW  CONDITIONS    IN   COMPRESSOR    STAGE   WITH   VELOCITY  DIAGRAM 
OF   FIG.  2  ON   CYLINDRICAL  STREAM  SURFACES 


DATA 


AT   M 

EAN   RADIUS   R       : 
in 

Am    * 

4-5° 

D    .   0. 

35 

X  =    . 

28  78  3     (E^.4) 
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.  41462^.  !•) 

E* 

*/«„ 
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\% 
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TABLE    31 

THREE-DIMENSIONAL   FLOW  CONDITIONS   IN   COMPRESSOR   STAGE  WITH  VELOCITY  DIAGRAM 
OF  FIG.    2  ON  CYLINDRICAL  STREAM   SURFACES 
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TABLE  12 

THREE-DIMENSIONAL    FLOW  CONDITIONS    IN   COMPRESSOR   STAGE  WITH   VELOCITY  DIAGRAM 
OF   FIG. 2   ON    CYLINDRICAL    STREAM   SURFACES 
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